Time evolution of the one-dimensional Jaynes-Cummings-Hubbard Hamiltonian 
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The Jaynes-Cummings-Hubbard (JCH) system describes a network of single-mode photonic cav- 
ities connected via evanescent coupling. Each cavity contains a single two level system which can 
be tuned in resonance with the cavity. Here we explore the behavior of single excitations (where an 
excitation can be either photonic or atomic) in the linear JCH system, which describes a coupled 
cavity waveguide. We use direct, analytic diagonalization of the Hamiltonian to study cases where 
inter-cavity coupling is either uniform or varies parabolically along the chain. Both excitations 
located in a single cavity, as well as one excitation as a Gaussian pulse spread over many cavities, 
are investigated as initial states. We predict unusual behavior of this system in the time domain, 
including slower than expected propagation of the excitation, and also splitting of the excitation 
into two distinct pulses, which travel at distinct speeds. In certain limits, we show that the JCH 
system mimics two Heisenberg spin chains. 



I. INTRODUCTION 

The realization of condensed matter models using non- 
traditional controllable systems has been an active topic 
of recent research. This is best exemplified by the work 
on cold atom quantum simulators [3, However the 
possibility of quantum optical simulators for condensed 
matter models has emerged. While finding links and 
analogies between different subfields of physics is inter- 
esting in its own right, of particular interest are effects 
that are difficult or impossible to observe in more con- 
ventional physical systems. 

There are a number of recent proposals [1, 0, H, 0, 
0, B B for constructing quantum optical condensed 
matter simulators. These use coupled cavity structures 
where confined photons are induced to interact via their 
coupling to embedded two-state systems. Possible two- 
state systems include color centers [1] , quantum dots [ill, 
[l3 | , superconducting strip- line resonators [H, [13, [3 , or 
coupled Rubidium microcavities (16] . 

In such systems, it should be possible to observe many- 
body effects such as quantum phase transitions where 
the particles of interest are photons (rather than elec- 
trons) and therefore can be readily injected, confined and 
observed individually. While many-body effects in the 
thermodynamic limit require a large number of cavities 
and/or photons (often solved in the mean field limit), it 
is quite clear that early experiments will be limited in 
both overall structure size and system controllability. It 
is therefore of interest to determine what (if any) effects 
can be observed in the few cavity/few excitation limit 
(where an excitation can either be photonic or atomic) 
[l9| . Importantly, the complex interactions 
between atomic and photonic components produce a va- 
riety of effects. Solitonic behavior has been predicted 
in the ID coupled cavity wave guide realizations of the 
Dickc model [20] and the XXZ model in the presence of a 
tilted magnetic field ^l| . If one is considering, for exam- 
ple, the type of localized/delocalized behavior typically 



found in strongly interacting systems, it is of vital im- 
portance not to be distracted by the background effects 
which come from the delocalized behaviour of photons 
themselves. For this reason, a detailed understanding of 
the single-particle dynamics is paramount. 

Here we investigate the time evolution of a one di- 
mensional coupled cavity waveguide [l^ described by the 
Jaynes-Cummings-Hubbard (JCH) model In the sin- 
gle excitation subspace we diagonalize the JCH Hamilto- 
nian exactly and consider the dependence of this system 
on three parameters: atom-cavity detuning, atom-photon 
coupling and coupling between cavities. We then investi- 
gate limits in which localized and delocalized behaviour 
can be seen. In this paper, we refer to the propagation 
of photonic or atomic excitations as the photonic and 
atomic components, or modes, respectively (in the atomic 
case, the atoms remain stationary whilst the excitation 
passes through). In particular we focus on three limits. 
Firstly, the limit where the atom-cavity detuning is zero 
and the coupling strength between cavities is much less 
than the photonic cavity-atomic coupling. In this limit, 
the propagation dynamics of the atomic and photonic 
components are identical; we find they propagate as a 
pulse travelling back and forth along the line of cavities. 
Secondly, we study the limit where photon hopping dom- 
inates all other parameters of the system. In this limit, 
the atomic component does not move, while the pho- 
tonic component propagates. Thirdly, we study the limit 
where the atom-cavity detuning is much larger than all 
other energy scales of the system. In this limit the atomic 
and photonic modes travel at two different speeds. In 
these three limits we compare the behaviour with that of 
two uncoupled Heisenberg spin chains [^, [H, [13] , in the 
one-excitation case. 

We discuss the JCH model and the uniform coupling 
case in Sec. [Hi We then continue with dispersion-free 
pulses using parabolic coupling, Sec. lIIIl and initial Gaus- 
sian pulses. Sec. |TVl In Sec. |V] we discuss the limit of 
large atom-cavity detuning, in this limit the atomic and 
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FIG. 1: (Color online) Visualizations of a ID JCH system. 
Fig. (a) shows an example of a coupled cavity waveguide in a 
photonic crystal. A lattice of holes in the membrane provide 
variations in refractive index, hence trapping photons of a 
given frequency. The red spheres indicate two-level systems, 
these are placed at sites without holes, called defects, which 
form photonic cavities. Fig. (b) is a visualization of a portion 
of a JCH chain as logical elements. The bottom row of circles 
indicate the photonic cavities, linked together by the hopping 
rate k. The top row of circles indicate the atoms, linked back 
to each photonic cavity by the coupling strength /3 [defined in 
Eq. ©]. 

photonic modes travel at two distinct speeds. 

II. UNIFORM COUPLING 

The JCH Hamiltonian (7i'"-^^) describes a system of 
A'^ cavities linked via photon hopping under the tight- 
binding approximation 0, [^, [l^ . A possible arrangement 
and schematic showing the couplings is shown in Fig. ([1]). 
The Hamiltonian is 

N N 

^acH^^^.c_,^^^^^t,^^ (1) 
i=i i,i=i 

where k is the inter-cavity coupling, and A is the adja- 
cency matrix, which is defined according to the geome- 
try and boundary conditions of the system. Using the 
geometry as implied by Fig. (HI) (a), that is, one spatial 
dimension with hard wall boundary conditions, A is given 
by 

1 0, otherwise ' ^ ' 

The Jaynes-Cummings Hamiltonian 7i^^ [2^ describes 
a single atom-cavity system. It can be written 

Hf^ = e^a+a:r + uj,a\a^ + |3^{(J+a, + CTra|), (3) 



where a| (a^) and cr^ are the photonic and atomic 

raising (lowering) operators respectively, and h = 1. In 
cavity i, the energy of the atom is given by e^, the cavity 
resonance by cui, and the cavity-atom coupling by f3i. 

The single cavity bare basis consists of states of the 
form |s, n) , where s £ {5, e} represents the two- level atom 
in the ground or excited state, and n is a non-negative 
integer, representing the number of photons in the cavity. 
The N cavity bare basis simply consists of the tensor 
product of N such single cavity bare bases. 

The single cavity dressed basis consists of the state 
\g, 0) (this state is a member of both the bare and dressed 
bases) and states of the form |±,n), where n is the 
number of excitations (photonic or atomic) in the cav- 
ity. They are energy eigenstates of the Jaynes-Cummings 
Hamiltonian, and are related to the single cavity bare ba- 
sis states by 

^ /3v^|g,n) + KA/2)±xH]|e,n-l) 

yVRrxRA (4) 
Vn > 1, 

where A = w — e is the detuning, and we have used the 
generalized Rabi frequency 

Xin) = \/n/32 + A2 /4 Vn > 1. (5) 

Consider a restricted basis of the full N cavity bare 
basis, consisting only of basis states that contain one ex- 
citation (atomic or photonic) in the whole system: the 
one excitation subspace. The total number of excitations 
in the JCH system is conserved, hence this is a restric- 
tion, not an approximation. We consider the case where 
all cavities are equal, hence Ci = e,uji = uj, /3i = f3 for 
all i = 1, . . . ,N. We write the restricted one excitation 
basis as \Q) (g) \g, 1) {\Q) |e, 0)), where Q e {1, . . . , N}, 
to mean a photonic excitation \g, 1) (atomic excitation 
|e,0)) at cavity number, Q, and \g,0) at every other cav- 
ity position. The states \Q) form a valid Hilbert space. 
In this restricted basis and under these conditions, the 
Hamiltonian ([1]) is greatly simplified, and can be repre- 
sented 

ni^^,^ JIN'S) z + PIN'S) x-KA(g>^-^^, (6) 

where /,„ is the m x m identity matrix, X and Z are 
the usual Pauli matrices acting on the atom-photon cav- 
ity subspace. The operator A in matrix form is given 
by the adjacency matrix of the connectivity graph, for 
example that given in Eq. In this form, the first 

and second terms affect only the atomic/photonic modes 
locally, but do not move the excitation to any other cav- 
ity. The first term describes the detuning and the second 
term the coupling between photonic and atomic excita- 
tion modes. The third term describes photonic coupling 
between cavities. 
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FIG. 2: (Color online) Space-time diagrams for evolution of an excitation (jl) ® {\g, 1) + je,0))/\/2) along a chain of 100 JC 
cavities for three different parameter regimes. In each case, we plot the probability of occupation of a particular cavity (vertical 
axis) as a function of time (horizontal axis) for the system initialized in an even superposition of the atomic and photonic mode 
of the first cavity. The upper three plots show the population of the atomic components whereas the lower three plots are the 
photonic components. The ratio of cavity-cavity coupling to atom-photon coupling varies from left to right, with (a) and (b) 
k//3 = lO""*, (c) and (d) = 10 and (e) and (f) = 10^. In plots (a), (b) and (f), the dashed lines represent Qa, as given 
in Eq. (0. 



By diagonalizing A the whole JCH Hamiltonian re- 
stricted to one excitation can be diagonalized. We now 
solve for the ID chain with hard wall boundary condi- 
tions. The eigenvectors of the A matrix given in Eq. ([2]) 
are fSl 



where the second term of this appears very similar to 
the Rabi frequency x(l)i Eq. ([5]), with the detuning A 
shifted by the cosine term. These eigenvalues have cor- 
responding eigenvectors 



|fc) 



V2(-l)^sin(^ 



lsin(^ 



N 
Q=l 



Qkn 
N+1 



\Q), (7) 



where k = 1, . . . ,N. 

We now wish to obtain the energy eigenstates and en- 
ergy eigenvalues for the entire one excitation subspace. 
For a linear chain of cavities the Hamiltonian Hj^^ 
can be expressed in the basis {\k) (g) \g,l), \k) |e,0)}, 
k ~ 1, . . . , N as a block diagonal matrix, in which the 
fcth block appears as 



The eigenvalues of the full Hamiltonian, from Eq. (O , 
are 



|±,fc) = 



(A + 2i;|)|e,0) +2/31.9,1) 
(A + 2^1)2 +4/32 



|fc). (10) 



This is an exact diagonalization of the linearly-coupled 
JCH system in the one excitation subspace. Hence, it is 
possible to determine analytically, for arbitrary N, the 
time-evolution of an arbitrary initial state. In the case 
with two cavities, i.e. N = 2, and in appropriate limits to 
each case, this evolution recovers equations (22), (28) and 
(29) of Ogden et al. Q. Specifically, we are interested in 
plotting the expectation value of both the photonic and 
atomic excitations in cavity j , by examining the number 
operators a^Oj and (^^^J ■ 

As a demonstration of the dynamics of the system 
when on resonance (i.e. A = 0), we consider evolution 
of an excitation initially located in the first cavity in a 
line of 100 cavities, in an equal superposition of atomic 
and photonic modes, i.e. 



Ej. ^ K cos 



kiT 

N+1 



± 



A f kn 



1 2 



(9) 



m^O)) - |l)®(|.g,l) 



,0))/%/2, 



(11) 



which corresponds to the single cavity energy eigenstate 
| + , 1) from Eq. ^ in the case A = 0. As such, there will 
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FIG. 3: Dispersion of the wave packet, AQ, at a fixed point in time, as a function of with atom-cavity detuning A = 
and number of cavities A*' = 100. The solid (dashed) line shows ^Qiu!^ic (AQphotonic)- The dotted (dot-dashed) line shows 
^qHcis ^j-j-jj J = 2k (J — k). The dispersion is measured at a time T = N/Ak, at which point the packet is 'freely' evolving 
along the chain. This point is chosen such that the effects of the hard-wall boundaries can be ignored for both large and small 
K. At the left of the plot, ^ 1, and the JCH chain mimics two identical Heisenberg spin chains, an example is shown of 
this localized behavior in Fig. ©(a) and (b), where k//3 = 10~^. The corresponding spatial profile of the pulse at T = N/4k, is 
shown in the left inset (the solid line shows the photonic profile and the dashed line shows the atomic profile, in this case they 
are coincident). At the right of the plot, 3> 1, and the photonic mode of the JCH chain mimics a Heisenberg spin chain, 
while the atomic mode does not propagate at all. An example of this type of localized behavior is shown in Fig. ©(e) and (f), 
where k/P = 10^, the corresponding profile of the pulse at T = N/Ak, is shown in the right inset. In the middle of the plot a 
travelling excitation shows a large amount of dispersion, an example of this delocalized behavior is shown in Fig. ©(c) and 
(d), where k//3 = 10, and the corresponding profile at T = N/4:K is shown in the middle inset. The horizontal lines show the 
dispersion AQ^^"^ of a Heisenberg spin chain with 100 spins after the initial state |1) evolves until the front of the excitation 
is half way (dotted line) and quarter way (dot-dashed line) along the excitation chain. The exact correspondence between the 
Heisenberg and JCH systems is seen in these asymptotic limits. 



be no evolution within a cavity, but only evolution be- 
tween cavities. We solve the evolution exactly, using the 
usual Schrodinger equation, idt\ip) = H\ip). This state 
also constitutes one of the simplest states to realize ex- 
perimentally, as the JC resonance can be driven directly 
by a transverse field. The dispersion and other effects dis- 
cussed in this paper stem directly from the strong atom- 
photon interaction via the JC Hamiltonian. If we were to 
consider a system with vanishing coupling to the atoms 
(/3 = 0), we regain the conventional photon propaga- 
tion results [13, [isj in which a Gaussian wave propagates 
smoothly. 

We are ultimately interested in the behavior of the sys- 
tem when various energies dominate, such as the atom- 
photon coupling or the cavity-cavity coupling. For this 
reason, we will now consider the system in several dif- 
ferent limits. When k//3 ^ 1 and A = 0, the atomic 
and photonic modes have identical propagation dynam- 
ics. The case k/(3 = 10~^ is shown in Fig. ©(a) and 
(b). We may understand the equal propagation because 
the atom-photon coupling /3 is much stronger than the 



cavity-cavity coupling k, hence the excitation is free to 
form the single cavity eigenstate between inter-cavity 
hops. By shifting to the interaction picture [13], that 



is Ix) 



the Hamiltonian Eq. ^ becomes 



H = — A<. 
2 



(12) 



where the fast rotating terms have been ignored. It is 
useful at this point to compare with the Hamiltonian that 
describes the well-known Heisenberg spin chain [2^, [2^ . 



N 



'1 



(13) 



n=l 
N r 



n=l 
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5Z QZ 



The single site basis of this system consists of spins 
pointing up and down along the direction of the z axis, 
{IT), li)}- The N site basis comprises of a tensor product 
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FIG. 4: (Color online) Evolution of a Heisenberg spin chain 
Hamiltonian [Eq. (IT4)l ]. with J = 1. (a) Evolution of |1) un- 
der uniform coupling, note that at first the excitation travels 
fairly neatly, but both within a single pass, as well as dur- 
ing the reflection, the excitation pulse spreads out. Note the 
faint lines which are parallel but displaced from the wavefront. 
These lines are due to the nature of the Heisenberg Hamilto- 
nian. (b) Evolution of |1) under parabolic coupling, note how 
the evolution is smooth and repetitive, the excitation trav- 
els sinusoidally from one end of the line of spins to the other 
end and back again, without spreading out. (c) Evolution 
of a Gaussian pulse as in Eq. (|26[) . with Qc ~ 50, s — 10, 
k — n/2, with uniform coupling. Note how the evolution is 
dispersion-free. 



of N such bases. This Hamiltonian also conserves the 
total spin in the z direction. So, if we limit the N site 
basis to having only one |t), and the rest ||), then the 
Hamiltonian can be represented in this restricted sub- 
space as (ignoring contributions from the term, which 
is largely a phase factor apart from a minor shift at the 
ends of the chain, this is inconsequential for long chains) 



H 



Hcif 



2 ' 



(14) 



where the basis vectors are now |1), |2), . . . , |A^), and \Q) 
represents |t) at site Q and |J,) at every other site, and 
A is the adjacency matrix, for example that which was 
introduced in Eq. The initial state |1) will evolve 
along the chain (due to the hard wall boundary condi- 
tions) with an approximate speed J (the rate at which 
the front of the excitation wave travels across the chain) , 
according the triangle wave 



N-l 



arcsm < sm 



Jt 



N+1 



(15) 

We find that the speed J as described above is given by 
the derivative of Eq. (fT5)) . i.e. J = \dtQ\ (where defined). 



FIG. 5: Evolution of a Heisenberg spin chain Hamiltonian 
[Eq. (I14|l ]. with J — 1. The solid lines display the expec- 
tation value of the position (Q), the dashed lines show the 
dispersion AQ. (a) Evolution of |1) under uniform coupling, 
AQ increases with time, this indicates increasing dispersion 
as time progresses, (b) Evolution of |1) under parabolic cou- 
pling, here AQ does not increase with time, as such this in- 
dicates dispersion-free evolution, (c) Evolution of a Gaussian 
pulse as in Eq. (l26|l . with Qc = 50, s = 10, k = 7r/2, with 
uniform coupling. As in (b), there is no overall increase in 
dispersion, as such this pulse is also dispersion-free. 



By comparing Eq. (|T4)) with Eq. (|12p . we see that the 
JCH Hamiltonian in the regime A = 0, k <C /3, mimics 
two Heisenberg spin chains, in each of the photonic and 
the atomic states. In both cases, J — k (the approximate 
excitation speed). 

Next we consider an alternative limit, where the cou- 
pling between cavities dominates the evolution, k//3 3> 1. 
In this regime, the atomic mode does not propagate at all, 
while the photonic mode propagates at twice the speed 
of the previous case. This is because the atom-photon 
coupling is effectively zero compared to the much faster 
inter-cavity coupling rate, freezing the atomic excitation. 
The case k/P = 10'^ is shown in Fig. ©(e) and (f). In 
this limit, the Hamiltonian trivially reduces to 



H 



-kA ( 



h + Z 



(16) 



as such, the photonic excitation travels with speed 2k, 
and the atomic excitation does not move at all. We can 
think of this limit as that of a pure photon gas, albeit 
with only one photon, which is equivalent to a single- 
excitation spin chain. 

When A = 0, and k and (3 are of the same order of 
magnitude, the evolution of the state Eq. p!T|) now expe- 
riences a large amount of dispersion. This is because the 
evolution no longer approximates two Heisenberg spin 
chains (as it does in the k//3 ^ 1 and k/P ^ 1 limits), 
but rather the full JCH nature of the evolution is ex- 
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FIG. 6: (Color online) Space-time diagrams for evolution of an excitation (|1) (g) {\g, 1) + |e,0))/\/2) along a chain of 100 JC 
cavities with parabolic inter-cavity coupling profile. The upper plots show the population of the atomic components whereas 
the lower plots are the photonic components. The ratio of cavity-cavity coupling to atom-photon coupling varies from left to 
right, with (a) and (b) k//3 = 10"", (c) and (d) k/P = 10°, (e) and (f) = 10^. In plots (a), (b) and (f), the white dashed 
lines represent {Q{t)}, as given in Eq. (|29p . The effect of the parabolic coupling is to constrain the pulse, resulting in well 
defined and reversible dispersion. 



pressed. As the atom-cavity coupling /3 and the cavity- 
cavity coupling K are of the same order, there is no oppor- 
tunity for the excitation to form any single cavity eigen- 
states, instead the excitation is relatively free to roam 
between photonic and atomic modes, as well as between 
cavities. The example k/P = 10 is shown in Fig. ©(c) 
and (d). 

To more clearly see this interplay between JC and pho- 
ton dominated regimes, we consider the dispersion of the 
wave packet as it travels along the chain. Fig. ^ shows 
how the dispersion of a pulse changes with We de- 
fine the position operators 

Qphotonic = diag(l, 0, 2, 0, . . . , iV, 0) 
Qatomic = diag(0,l,0,2,...,0,iV) ^ ' 

in the basis {|(5)(g)|g,l),|(3)«)|e,0)},Q = 1,...,7V. The 
definitions of Eq. ^Tl} assume that when calculating ex- 
pectation values for photonic (atomic) position, the pho- 
tonic (atomic) states are selected from the wave function 
and normalized as a single vector of length N [s^l • 

The dispersion for photonic and atomic modes is given 
by the standard deviation of the position operator AQ 

(AQtype)' = (Qt'ypc) " (Qtype)^ (18) 

Looking at Fig. ([3]), we see the link between the two 
limits k/P < 1 and k/(3 > 1 for A 0. The sohd 
(dashed) hne shows AQig^,^ (AQJ^o«to„ic)- In these two 
limits, the JCH chain mimics two Heisenberg spin chains. 



while in the middle range, the full JCH dynamics are re- 
alized. With increasing k//3, the photonic (atomic) dis- 
persion of a pulse is constant (constant) when n/P <^ 1, 
increases for moderate values of k//3, then becomes con- 
stant (zero) when k/(3 ^ 1. Dispersion is measured at 
the time T = N/Ak for a system with N = 100 cavities 
initially in the state given by Eq. (fTT|) (i.e. |1) (g) \+, 1) 
with A = 0). The horizontal lines drawn on the figure 
show the dispersion AQ^ois (gHois ^ diag(l, 2, . . . , TV)) 
of a Heisenberg spin chain with 100 spins and initial state 
|1) for J — K (dotted hne) and J — 2k (dot-dashed hne) 
along the excitation chain. The J — k Heisenberg dis- 
persion line matches both the photonic and atomic dis- 
persion in the limit when k/ f3 <^ 1. The J = 2k Heisen- 
berg dispersion line matches the photonic dispersion in 
the limit when k/(3 ^ 1, while the atomic dispersion in 
this limit tends to zero (as the atomic mode no longer 
propagates). 

The left inset of Fig. ^ and Fig. ©(a) and (b) corre- 
spond to the example k/ (3 ^ The inset shows the 
pulse profile at the time T — N/Ak [with pulse approxi- 
mately a quarter along the chain, = N/4, and sohd 
(dashed) line representing the photonic (atomic) profile] , 
while part (a) and (b) shows that the evolution of the 
atomic and photonic modes are identical. The middle 
inset of Fig. ([3]) and Fig. ©(c) and (d) correspond to 
the example = 10. The inset shows the high dis- 
persion of the pulse profile at the time T — N/Ak (the 
dashed line shows atomic profile), and part (b) shows the 
high dispersion evolution. The right inset of Fig. ^ and 
Fig. ©(e) and (f) correspond to the example k/ (} = 10^. 
The inset shows the pulse profile at T = N/4k, note 
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FIG. 7: Dispersion of the wave packet, AQ, at a fixed point in time for a JCH system with parabolic coupling. The solid 
(dashed) line shows AQf^mic (A(3p£<5J;onic)- The dotted (dot-dashed) line shows AQ^"'' with J = 2k {J = k). The parameter 
ranges are identical to Fig. ([3]), although the range of the vertical axis is different. The insets show the pulse profile at three 
values of The solid (dashed) line represents the photonic (atomic) profile. The result is that, while the dispersion is 

considerably larger than the uniform coupling case, in this case the pulse is more 'well behaved'. More precisely, the pulse has 
a Gaussian profile away from the boundaries and its evolution is completely reversible. In the extreme limits for k, we have 
again perfect agreement with that expected from the appropriate spin chain model. 



that the photonic profile is approximately half way along 
the chain (Qa = while the dashed line shows the 

atomic profile, which does not propagate. The atomic 
and photonic modes remain in superposition, despite the 
fact that the atomic mode does not propagate and the 
photonic mode does. 



III. PARABOLIC COUPLING 

One way to control the dispersion of the pulse is to 
define a non-uniform distribution of couplings between 
cavities. One such possibility is to choose a regime of 
parabolic couplings [2^, , where the coupling between 
cavity i and cavity i -I- 1 is given by \/i{N — i), so that 
the adjacency matrix becomes 



Consider the profile (probability distribution at an in- 
stant in time) of a single excitation as it travels along the 
chain. In Fig. (01) (a), we display the expectation value of 
the evolution of a single up spin (initially at |1)) for a 
Heisenberg spin chain with uniform coupling J between 
each of the 100 sites. One sees that the pulse is initially 
well-formed, but later suffers from increasing dispersion 
of the excitation, both when travelling through the chain, 
and also when reflecting from the end of the chain. Re- 
gardless, the speed of the pulse is given by J, at least for 
the first reflection, where the wave speed can be realisti- 
cally interpreted. Fig. ([51) (a) is a different way of show- 
ing the same evolution, by plotting the expectation value 
and uncertainty of the operator Q^""* = diag(l, 2, . . . , N) 
highlighting the spread of the pulse with time. While the 
physics of Heisenberg spin chains is well studied [H, [2^ , 
it is useful to contrast this behaviour with the next regime 
we will examine. 



Ai,, = { ^i[N -i) j -i = l 
otherwise 



(19) 



Thus the couplings will be symmetric around the cen- 
tral cavity (or between the two central cavities), with 
strongest coupling at the center of the chain and weakest 
coupling at the ends. A spin chain of length N with this 
coupling can be mapped to a single spin s — {N — l)/2 
particle, placed in a magnetic field in the x direction. 
This system provides the phy sical insight for why this 
coupling is dispersion- free [30|- The eigenvalues of this 
matrix are Ek = N — 1 — 2k and the eigenvectors are 



\k) 



N 

E 

Q=l 



il-N)k{N-l)\ 



(-l)'=fc!2^-i((3- l)!(iV-Q)! 



KkiQ-l,-,N-l)\Q), 



(20) 
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where the Pochhammer symbol {N)k is defined as 



for a Heisenberg spin chain, or 



(iV)fc = iN)iN + l)...{N + k-l), k = l,2,3,... (21) 

and the Krawtchouk polynomial K [3l|, [12] is related to 
the hypergeometric function F by 



Kk{l,P,N) = 2F1 (■ 



-k,-l\ -1) 



(22) 



Now we examine the behaviour of the JCH chain un- 
der this parabolic coupling scheme. We find that, again 
in the basis given by {|fc) (g) \g, 1), \k) (g) |e, 0)}, the Hamil- 
tonian is a block diagonal matrix, with the fcth block 
appearing as 



k{N - 1 



2k) (i 
-A/2 



The eigenvalues of the full Hamiltonian 
parabolic coupling case are 



(23) 



in the 



]^[n{2k 



N + 1) 



± V[A + K(2fc -N + 1)]2 + 4/32 1 (24) 
and the eigenvectors are 



|±,fc) 



(A + 2g^)|.g,l)+2/3|e,0) 
(A + 2£:|)2 + 4/32 



® \k). (25) 



We find the matching that occurs in limiting cases for 
the parabolic JCH system, occurs in exactly the same 
way as the uniform coupling case. That is, when A = 
and K <^ (3, J = K for both the atomic and photonic 
parts. When A = and k:$> (3, J = 2k for the photonic 
mode and the atomic mode does not propagate. Figs. ^ 
and ([7]) show very similar results to Figs. ^ and with 
only one significant difference (apart from the dispersion- 
free pulses), that the uncertainty in position for moderate 
values of k/ f3 is much larger than the uniform coupling 
case. 



IV. GAUSSIAN PULSES 

As an alternative to modifying the coupling profile, 
we can consider uniform coupling and a Gaussian wave 
packet as our initial state. In this case, the momen- 
tum distribution of the pulse is well defined (and nar- 
row) which allows the excitation to travel down the chain 
with minimal increase in dispersion [13, [3 • We therefore 
choose an appropriate initial state 



N 



IV' 



(fc,S,Qc) 



-ikQ 



\Q). (26) 



|^(fc,.,Qe)(^ = 0)) 

N _ 2 

= ^ e-^^^e-'^^IQ) ® (|g, 1) + |e, 0))/^/2, 
Q=l 

(27) 

for a JCH chain, where M is the normalization, k is the 
wave number, Q,. is the center of the pulse, s is the width 
of the pulse, and Q denotes cavity number, with Q — 
1, . . . ,N. In this paper we choose Qc — N/2, s — N/10, 
such that the pulse is initiated at the center of the chain 
with a width approximately 2/10 times the length of the 
chain. The value k ~ tt/2 + wk, where n is an integer, 
produces dispersion-free evolution of the pulse. 

Fig. ([H]) shows the evolution of a Gaussian pulse for 
three different values of k//3, in analogy to Fig. ([2]) and 
^ . The white dashed lines correspond to a triangle wave 
similar to Eq. (fT5|) . but with a phase shift by 7r/2. The 
phase shift is necessary as we have initiated the Gaus- 
sian pulse in the center of the chain, to avoid boundary 
effects (contrary to the first uniform coupling case, where 
the boundary effects were essential for the motion of the 
pulse). The dashed fine is given by 



Q Gaussian 
A 



N -I 



■ arcsm 



■ . f J-Kt 

sm 

. V N 



N +1 



(28) 



In the limits k//3 <C 1 and k//3 ^ 1, the Gaussian pulse 
with the choice of k — 7r/2 evolves under simple transla- 
tion motion, interrupted by reflection off the boundaries. 
This translational motion, as well as the interference pat- 
tern as it reflects off the boundary, mimics exactly the 
propagation of a Gaussian wave packet under the one- 
dimensional Schrodinger equation [13, [H, HI] . 

Fig. © shows the dispersion of the pulse as a function 
of k//3, in analogy to Fig. ^ and ([7]). Note how this fig- 
ure is much simpler than the earlier two, as the dispersion 
is a constant throughout the evolution of the pulse, pro- 
vided the pulse is not interacting with the boundaries. 
As such, only one horizontal line to indicate dispersion 
of a corresponding Heisenberg spin chain is necessary. 
Also, the dispersion of the photonic mode is equal to the 
dispersion of the atomic mode. 

By comparing Fig. (HI) (a) (uniform coupling) with 
(|4])(b) (parabolic couphng) and Fig. ^{c) (uniform cou- 
pling, Gaussian pulse) we can indeed see that the pulse 
in the uniform chain has more dispersion than that of 
the parabolic chain and the uniform chain with an initial 
Gaussian pulse, the dispersion is manifest as a number 
of faint lines parallel to the main wavefront. Further, in 
Fig. ([5]) we have plotted the expectation value of the po- 
sition {Q) (solid lines) and AQ (dashed lines), for both 
the uniform coupling (a) , parabolic coupling (b) and uni- 
form coupling with Gaussian pulse (c) cases. In fact, the 
expectation value of the position Q in the parabolic chain 
is of this simple form: 



9 



Expectation value of atomic and photonic excitations 




100 200 300 100 200 300 100 200 300 400 

Time (in units of k'^) 



FIG. 8: (Color online) Space-time diagrams for evolution of an Gaussian pulse along a chain of 100 JC cavities with a uniform 
inter-cavity coupling profile. The upper plots show the population of the atomic components whereas the lower plots are the 
photonic components. The ratio of cavity-cavity coupling to atom-photon coupling varies from left to right, with (a) and (b) 

= 1Q"^ (c) and (d) = 10°, (e) and (f) k/(5 = 10^ In plots (a), (b) and (f), the white dashed lines represent 
as given in Eq. (|28|l . The effect of the Gaussian pulse is to have chosen a narrow momentum distribution, which allows the 
excitation to propagate freely in the system. 



{Q[t)) = -[N +1- [N -I) cos Jt]. (29) 

Hence the period of oscillation is 27r/J, which does not 
depend on the number of cavities N. For the cou- 
pling profile, we find that the pulse (while more dis- 
persed) is approximately Gaussian shaped in space. This 
can be seen by examining Fig. (fTUl) (b). which shows 
the pulse profile at fixed instants in time. In com- 
parison. Fig. (fTOl) (a) has much higher dispersion, while 
Fig. (fTm) (c) also has low dispersion, which is largely con- 
stant and only reduces slightly as the pulse is reflected 
at the boundaries of the chain. 



V. LARGE DETUNING LIMIT 

We now consider the case where the magnitude of the 
detuning is much larger than the other energy scales of 
the system, that is, k, /3 <C |A|. The effect of increas- 
ing detuning is to decrease the atom-photon coupling, 
thereby approximating a Bose-Hubbard system. We are 
therefore interested in any non-trivial effects which ap- 
pear in this limit. 

We now shift again to an interaction picture, where 
states of the system are transformed to 

Under this shift, the Hamiltonian becomes 



, \i + z f 2/32 \ 
L 2 V A2 + 4/32; 

I-Z 2/3^ A/3 1 

2 A2 + 4/32 A2 + 4/32 .■ 

In this equation, the coefficient of (/ + Z)/2 is related 
to the speed of the atomic component, and the coefficient 
of {I — Z)/2 is related to the speed of the photonic com- 
ponent. We will discuss the effect of the X term shortly. 

To explore the dynamics of this system, we again 
start with the equal superposition initial state as given 
in Eq. For example, the uniform coupling case 

with A = 10^/3, K = /3 is shown in Fig. CI]) (a). By 
again comparing to the Heisenberg spin chain Hamil- 
tonian, Eq. (fT4)) . we see now that the photonic excita- 
tion corresponds to a Heisenberg spin chain with speed 
J = 2k[1 - 2/32/(4/32 + A2)], and the atomic excita- 
tion corresponds to a Heisenberg spin chain with speed 
J — 2k/32/(4/32 -I- A2). In this regime, the atomic and 
photonic excitations travel at two completely different 
speeds: the photonic mode travels much faster than the 
atomic mode, although the integrity of each mode is pre- 
served. We also can see in Fig. (|TT|) that this splitting 
of the system into two separate modes happens for both 
parabolic coupling (b) and an initial Gaussian pulse (c). 
This splitting is an observation which can be predicted di- 
rectly from the kA independent form of the various com- 
ponents of Eq. (ini). Note that the final term of Eq. (pij) 
is significantly different than both the time scale of the 
photonic propagation (w 2k) and the time scale of the 
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Uncertainty in position at T=N/8k (uniform coupling, Gaussian pulse) 




-3-2-1 1 2 3 4 

iogio(K/p) 

FIG. 9: Dispersion of a Gaussian wave packet, AQ, at a fixed point in time, as a function of hi/ P, at A = 0. Tlie dispersion is 
measured at time T = N/8k, at wliicli point the packet is 'freely' evolving along the chain. The initial state of the Gaussian 
is given by Eq. (|27p . with center of the pulse at Qc ~ N/2, width of the pulse s = A'^/10, wave number k — 7r/2, and number 
of cavities A'^ = 100. The width of the pulse is small enough such that the boundaries do not have an effect on AQ for the 
value of T = N/8k, in the limits k;//3 <^ 1 and ^ 1. At the left of the plot, /t//3 <^ 1, and the JCH chain mimics two 
identical Heisenberg spin chains, an example is shown of this localized behavior in Fig. ©(a) and (b), where k//3 = 10~^. The 
corresponding spatial profile of the pulse is shown in the left inset, (the solid line shows the photonic profile and the dashed 
line shows the atomic profile, in this case they are coincident). At the right of the plot k/ P 2> 1, and the photonic mode of 
the JCH chain mimics a Heisenberg spin chain, while the atomic mode does not propagate at all. An example of this type 
of localized behavior is shown in Fig. (I8])(e) and (f), where k//3 = 10'^, the corresponding profile of the pulse at T = N/8k is 
shown in the right inset. In the middle of the plot a travelling excitation shows a large amount of dispersion, an example of 
this delocalized behavior is shown in Fig. (|8}(c) and (d), where k/P = 10, and the corresponding profile at T = N/Sk is shown 
in the middle inset. The horizontal line shows the dispersion AQ^^"^ of a freely evolving pulse [this value does not change, 
provided the pulse is not interacting with the boundary, as is evidenced by Fig. (I5])(c)]. 



atomic propagation (w 2k/5^/A^). As such it does not 
contribute significantly to either mode of propagation in 
the very large detuning limit. 

If the detuning is reduced such that this separation of 
time scales is not so strong, we must ask the question 
what is the effect of this final X term in Eq. ([3T|) ? To 
answer this, we write it out in the original basis of atoms 
and photons on neighboring sites j and j + 1, 

^j.j+i = + 0-74+1 + '^i+i^J + ^i+A^ (32) 

in which we see that it is an effective JC type coupling 
between neighboring atoms and photons. The effect of 
this term on the propagating atom and photon chains 
will be zero at first order as it couples between chains. 
The second order effect is not zero, as 

/ + - + - I t t ^ ^ > 

giving an effective next-nearest-neighbor coupling (as- 
suming the one excitation subspace). The strength of 



this correction is given by the square of the X coef- 
ficient in Eq. (|3ip and therefore vanishes as the detuning 
is increased. In the limit where mixing effects are visible 
due to this additional term, we also observe asymmetric 
behavior with respect to the sign of the detuning, simi- 
lar to that seen previously in a two cavity system Q. It 
should be noted that here, we always initialize the sys- 
tem in an equal superposition of atomic and photonic 
states, whereas in previous work, the detuning asymme- 
try is accentuated by an initial state which is always an 
eigenstate of the JC Hamiltonian. 



VI. CONCLUSIONS 

In conclusion, we find that when limited to the one 
excitation subspace, one can analytically solve for the 
evolution of the excitation pulse. As such, we extend 
previous work on the topic [1, Q from TV = 2 to the 
many-cavity regime. We observe both localized and de- 
localized behavior in this system, which points to a com- 
plex interplay between atomic and photonic degrees of 
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(a) Pulse profile at various points in time (uniform coupling) 

' ' ' 1 

t=0 4t=N/4J 

.t=N/2J t=3N/4Jt=N/>^ 



Expectation value of atomic excitation state 
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Cavity number 

(b) Pulse profile at various points in time (parabolic coupling) 

t=0 t=7t/Ji 

=7i/4J , ,^ , t=3n/4J 
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(c) Pulse profile at various points in time (Gaussian pulse) 




40 60 
Cavity number 

FIG. 10: The pulse profile for various points in time as a 
function of cavity number, for evolution of a Heisenberg spin 
chain, (a) Evolution of jl) with uniform coupling between 
spins, (b) Evolution of |1) with parabolic coupling between 
spins, (c) Evolution of a Gaussian pulse, as given in Eq. 
with uniform coupling between spins. In the first uniform 
case, the profile is initially given by a Kronecker delta func- 
tion, but later spreads into a function with one primary peak 
and a number of smaller trailing peaks. In the parabolic case, 
at each end of the chain, the pulse is given by a Kronecker 
delta function, while in the middle of the chain the pulse ap- 
proximates a Gaussian. In the uniform coupling case with an 
initial Gaussian pulse, it evolves along the chain with fixed 
profile, and at the ends of the chain interferes and changes 
direction. 



freedom, even in the single excitation limit. 

We consider three natural limits, and show that the 
behaviour of the atomic and photonic modes of the JCH 
chain can be mapped to two independent Heisenberg spin 
chains. In the limit when the detuning is zero and the in- 
tercavity coupling is much smaller than the atom-photon 
coupling, we find that the system is mapped to two 
Heisenberg spin chains both with J = k. When the de- 
tuning is zero but atom-photon coupling is much smaller 
than intercavity coupling, the photonic mode propagates 
with approximate speed J — 2k and the atomic mode 
does not. 

In the limit when the detuning is much larger than 
both the intercavity coupling and the atom-photon cou- 
pling, the system is mapped to two Heisenberg spin 
chains, with J = 2k[1-2kI3'^ /{A'^+Af]'^)] for the photonic 
mode, and J — 2k/3^/(A^ -j- 4/3^) for the atomic mode. 
We also derive similar analytic solutions for the case of a 
parabolic variation in inter-cavity coupling, resulting in 
a Gaussian like wave packet propagation along the chain, 
as well as an initial Gaussian pulse in a uniform coupled 
chain. 

The JCH system provides an interesting playground 
for studying many-body physics in an atom photon con- 
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FIG. 11: (Color online) These plots show the evolution of 
the state |1) ® (|e,0) -I- |(?, 1)) (a) and (b) and a Gaussian 
pulse (c) in a one-dimensional JCH system consisting of 100 
cavities, with k = f3 and A//3 = 10''. Note the time scales 
are different in each case, (a) Evolution in a uniform chain, 
where the dashed lines are a triangle wave giving approximate 
evolution of the wavefront, as given in Eq. (llSp . (b) Evolution 
in a parabolically coupled chain, where the dashed lines are 
as given in Eq. (|29|l . (c) Evolution in a uniformly coupled 
chain, with an initial Gaussian pulse, as given in Eq. (|27|l . 
The dashed lines show the expectation value of position, as 
given in Eq. ([28|) . 
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text. Given that initial experiments will be limited in 
both cavity and excitation number, it is important we 
understand these temporal dynamics in a range of sys- 
tem size regimes. 
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